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both standards. 
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1. Introduction 

Optical atomic frequency standards based on single trapped and laser-cooled ions have 
made important progress and may soon reach and surpass the level of accuracy that is 
achievable with cesium atomic clocks in the microwave region today El El EI El El E| • 
An ion trap provides a very well-controllable environment and allows one to work at 
the conceptually simplest level: with one single atom. It is this simplicity that provides 
the basis for the high accuracy that can be obtained in spectroscopic measurements 
on these systems. 

In atomic frequency standards the frequency of an oscillator producing an 
electromagnetic wave is stabilized to an atomic reference transition. The error signal 
for the frequency lock is derived by modulating the oscillator frequency around 
the atomic resonance and by measuring the resulting modulation of the frequency- 
dependent excitation probability p to the upper atomic level. If only a single ion is 
used, a difficulty arises from the weakness of this signal and the consequently limited 
signal-to-noise ratio. The ion is probed repeatedly by pulses of the radiation. The 
state of the ion after an excitation attempt can be determined by applying laser 
radiation to induce resonance fluorescence on a transition that shares the lower state 
with the reference transition. This scheme was proposed by Dehmelt and is sometimes 
called electron shelving [HJ. Let us denote the two levels that are connected by the 
reference transition as \A) and \B) and assume that the ion is initially in the lower 
state | A). After an excitation attempt the ion generally will be in a superposition state 
a | A) + P\B) and the measurement with the electron shelving scheme is equivalent to 
determining the eigenvalue P of the projection operator P = \B) (B\. If no fluorescence 
is observed (the probability for this outcome being p = \(3\ 2 ) the previous excitation 
attempt is regarded successful (P — 1), whereas the observation of fluorescence 
indicates that the excited state was not populated (P = 0). In one measurement 
cycle only one binary unit of spectroscopic information is obtained. Under conditions 
where the average excitation probability p is 0.5, the result of a sequence of cycles 
is a random sequence of zeros and ones and the uncertainty in a prediction on the 
outcome of the next cycle is always maximal. These population fluctuations and 
their relevance in atomic frequency standards were first discussed by Itano et al., who 
named the phenomenon quantum projection noise (QPN) [5]. A simple calculation 
shows that the variance of the projection operator is given by |5] 



If the signal-to-noise ratio of a state measurement with a single atom is estimated as 
SNR = pj \Jp{\ — p), it can be seen that for p = 0.5, SNR = 1 only. 

The stability and the accuracy of an atomic frequency standard are ultimately 
determined by the signal-to-noise ratio, by the linewidth Ai/ and by the fluctuations of 
the center frequency vq of the atomic resonance The instability that is obtained 
after an averaging time T is conveniently expressed by the Allan deviation cr y (T) which 
can be written as 



This expression describes the case of a standard that is operated in interrogation 
cycles of duration t c . SNR denotes the signal-to- noise ratio that is obtained in one 
cycle and C is a numerical constant of order unity. The instability decreases inversely 
proportional to the square root of the number of cycles, giving rise to <J y (T) oc 1/y/T. 



(AP) 2 =p(l-p). 



(1) 




(2) 
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Parameters that are obtainable with present single-ion optical frequency standards 
@] 03 E] are of the order of magnitude Ai//^ = 10~ 15 and t c = 1 s. There are 
prospects to reach an accuracy in the range of 10 -18 |B1I2|, but from equation (2) it 
can be seen that with C/SNR « 1, an averaging time of about 10 6 s or roughly 10 
days would be required to complete a single frequency measurement at this accuracy. 
For the practical use of such a standard it will therefore be of high importance to 
find the interrogation conditions and servo techniques which minimize C/SNR, i.e. 
which make the most efficient use of the single-ion signal in order to provide a stable 
reference frequency. 

In this paper we present a theoretical and experimental study of the stability of 
a single-ion frequency standard under conditions where QPN is the dominant noise 
source. Because QPN presents a fundamental noise limit if only a single ion is probed, 
it defines the target below which technical noise due to the short-term frequency 
instability of the probe laser should be reduced. In the first part of the paper we 
present results of numerical calculations that allow one to find optimized excitation 
parameters for the case that the linewidth of the laser is smaller than the natural 
linewidth of the atomic reference transition. The choice of this regime is motivated 
by the fact that most optical frequency standards presently use reference transitions 
with natural linewidths in the range of a few hertz to a few hundred hertz pQ , whereas 
the most advanced laser sources already have subhertz linewidths ^J. With the 
practically unlimited interaction time for the trapped ion, this situation is in contrast 
to atomic microwave frequency standards where the natural linewidth is usually much 
smaller than the linewidth of the interrogation oscillator and where the resolution is 
limited by the finite interaction time. We compare the instability that is obtainable 
by excitation with a single laser pulse (Rabi excitation) and with two pulses (Ramsey 
excitation). The influence of a non-zero laser linewidth and of cycle dead time that 
is needed for the preparation of the ion and for state detection is discussed. While a 
final optimization will also have to take the specific frequency noise properties of the 
probe laser into account, these general results are meant to provide initial guidance 
in the selection of suitable experimental parameters. During the preparation of the 
manuscript two related studies were published |121 I13| and we will comment on the 
relation to that work at the end of section (3). 

In the second part of the paper we present a servo algorithm that can be used 
to stabilize the laser frequency to the single-ion signal. We investigate the dynamic 
behaviour of the servo system and the resulting stability of the frequency standard 
by numerical Monte Carlo simulations. Possible servo errors due to laser frequency 
drift are also considered. Results from the simulations are compared to data from 
experiments with two Yb + frequency standards that implement this servo algorithm. 
Previously, brief descriptions of frequency stabilizations to a single ion were given 
in El and a more detailed discussion including analytical expressions for the 
expected instability was presented by Bernard et al. |16| . 

2. Comparison of Rabi and Ramsey excitation 

Using numerical calculations, we study possible excitation conditions for the 
interrogation of the single-ion resonance in order to find parameters that yield 
a minimum instability. We make the assumption that the obtainable resonance 
linewidth and excitation probability are limited by the finite lifetime r of the excited 
atomic state. This regime will be typical for the operation of many single-ion optical 
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frequency standards that are presently studied |SJ 1171 I18| . In optimizing the 
excitation conditions, a compromise has to be found between the aims of obtaining 
a narrow linewidth and a high excitation probability. Because of their different 
dependence of the signal on the probe time i, we compare the two methods of Rabi 
excitation and Ramsey excitation. 

In order to calculate the shape of the excitation spectrum we use the vector 
representation of the optical Bloch equations in the form: 

u = Auv- i(r + r L )u (3) 

v = — Auju — uirw — — (r + Tl)v (4) 

w = ujrv - T(w - -) (5) 

where 2w — 2p — 1 is the population difference between the excited and the ground 
state and u and v are the real and the imaginary parts of the coherence, ljr denotes 
the Rabi frequency. The detuning of the laser frequency / from the atomic resonance 
vq determines Alo = 2ir(f — Vq). The damping rate T = 1/t describes the relaxation 
through spontaneous atomic decay and the additional coherence relaxation rate is 
used to model the influence of the laser linewidth. Tl describes a loss of coherence 
due to white laser frequency noise which would result in a Lorentzian excitation 
spectrum of full width at half maximum (FWHM) Tl/2it for Rabi excitation under 
the conditions u R < T L , T < T L , and 1/t < T L . 

By assigning a fixed value of the Rabi frequency ujr to the optical transition we 
neglect the influence of the finite temperature of the ion. This is justified as long 
as the product of the ion's mean vibrational quantum number (n) and the squared 
Lamb-Dicke parameter is small: (n)k 2 h/2mfl -C 1 (k: laser wavenumber, m: mass, 
Q: vibration frequency of the ion). This condition can be fulfilled by the use of a trap 
with high vibrational frequency or by cooling close to the vibrational ground state 
(n) ~ 0. Otherwise, the n-dependence of the Rabi frequency will wash out population 
oscillations and will make it impossible to deterministically excite the ion to state \B) 
by the application of a resonant 7r-pulse |12j . 



2.1. Rabi excitation with a pulse of duration t 

For various values of the Rabi frequency ujr and the pulse duration t, the Bloch 
equations were integrated numerically to calculate the maximum excitation probability 
Pmax at resonance. The width (FWHM) of the line profile Aojq was determined by 
searching for the detunings where the excitation probability drops to p max /2. We 
assume that these detunings are selected for the interrogation of the atomic line in the 
operation of the frequency standard. The suitability of the excitation conditions can 
then be determined by calculating a dimensionless stability parameter that we define 
as 



q 1 n 1 ^ Aujqt t + td 

dP = \ 7TPmax(l - -p max ) \/ . (6) 

V 2 2 p max V t 

With this definition, the Allan deviation (cf. equation (2)) is approximately given by: 

^- Sp 2kS- ^ 
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The first factor in Sp describes QPN for the interrogation at the half-linewidth 
detuning (cf. equation (1)) and the second factor is inversely proportional to the 
slope of the obtainable error signal. The last factor describes the expenditure of time 
for one excitation cycle, t c — t + td, consisting of the probe time t and the dead time td 
that is needed for detecting the state of the ion, for laser cooling and the preparation 
of the next interrogation cycle. The atomic lifetime r serves as the scaling parameter 
for time and frequency. The numerical search for the minimum of Sp was performed 
on a grid of values for t and lor that was chosen to determine the optimized parameters 
with a relative uncertainty of 2%. 

2.2. Ramsey excitation with two short pulses separated by a time t 

In the optimization of Ramsey excitation with two short pulses (of a duration that 
is much smaller than t), the Rabi frequency does not have to be treated as a free 
parameter because, irrespective of the relaxation rates T and Tl, the choice of two 
7r/2-pulses always ensures the highest signal contrast and lowest instability. In analogy 
to equation (6), the stability parameter for the case of Ramsey excitation is calculated 
as 



Here, p ma x and p m i n are the maximum and minimum excitation probabilities near the 
central Ramsey fringe and p = (p max +p m in)/2 is the average excitation probability. 
The full fringe width Aw is determined numerically as the frequency difference 
between resonance and the first minimum. 

By replacing Sp by Sp in equation (7), the Allan deviation obtainable with 
Ramsey excitation can be calculated. Since the actual slope of the spectral profiles 
at the FWHM points enters Sp and Sp only in the approximation Ap/ Aluq, we can 
not expect to obtain exactly identical Allan deviations in experimental realizations of 
the Ramsey and Rabi schemes with Sp = Sr, but we expect the two systems to show 
comparable instabilities to within a few percent under this condition. 

3. Optimized parameters for Rabi and Ramsey excitation 

Using the stability parameters defined in equations (6,8) as figures of merit, figure 
(1) shows the comparison of the optimized signal contrasts, linewidths, and stability 
parameters for Rabi and Ramsey excitation. It is assumed that is negligible 
compared to the natural decay rate T and that the cycle dead time td is much smaller 
than the probe time t. If the duration of the probe time is increased, a narrower 
resonance line is obtained but with a weaker signal. In the case of Ramsey excitation 
with two 7r/2-pulses, the signal contrast (figure (la)) drops exponentially like 



while the resonance linewidth (figure (lb)) narrows continuously like 1 ft and may reach 
values well below the natural linewidth Avq — 1/(2ttt). For Rabi excitation with an 
optimized pulse area close to the 7r-pulse condition ujpt — it, the excitation probability 
decreases slower than exponentially (in the limit of weak continuous excitation one 
would expect p m ax oc u) R ) but the linewidth never becomes smaller than A^o- The 
optimization of the stability for each value of t leads to pulse areas that are slightly 




(8) 



e -t/2r 



(9) 
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Figure 1. Comparison of the obtainable stability with Rabi (solid line) and 
Ramsey (dot-dashed line) excitation for the case of negligible laser linewidth. 
The optimized signal contrasts p ma x — Pmin (where p m i n = for Rabi excitation) 
(a), linewidths (b) and stability parameters (c) are plotted as a function of the 
probe time, measured in units of the atomic lifetime t. 



larger than for a 7r-pulse. In the range 0.5 < t/r < 5 the optimal pulse area can be 
approximately described by URt/n — 1.071 + 0.068 t/r. 

For Rabi excitation the lowest value of the stability parameter is Sp — 3.22, 
obtained for t — 1.88r and ojr — 2.02I\ For Ramsey excitation Sr = 2.59 can be 
achieved at t = I.OOt, indicating an advantage of about 20% for the optimized stability 
in comparison to the Rabi method. The excitation spectra for these two cases are 
plotted in figure (2) , showing that the optimization leads to spectra with quite similar 
shapes of the central resonance. Using equation (7) these results indicate that for the 
parameters of the Yb + frequency standard based on the 2 S\/2 — > 2 D^/ 2 transition 
with v a = 6.88- 10 14 Hz and r = 0.051 s, an Allan deviation o- y {T) = 2.7 • l^ 15 / ^/T/l 
should be achievable with Ramsey excitation. 

If realistic experimental conditions like non-zero laser linewidth and cycle dead 
time are taken into account, the stability difference between the two excitation 
methods becomes less pronounced. For the Ramsey scheme the additional relaxation 
Tl due to the laser linewidth does not broaden the resonance but reduces the fringe 
contrast. In the case of Rabi excitation the line is broadened and the maximum 
excitation probability is reduced. In both cases it is advantageous to reduce the probe 
time t below the values found for T L <C 1/r. For a laser linewidth that is equal to 
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Figure 2. Excitation spectra that yield optimal stability in the case of Rabi (solid 
line) and Ramsey excitation (dot-dashed line), corresponding to the minima of S 
in figure (lc). 



the natural linewidth (Tl/27t = A//n) the optimized parameters shift to t — l.lOr, 
ujr = 3.27r, and Sp = 4.28 for Rabi excitation, and to t = 0.50r and Sr = 3.64 for 
Ramsey excitation, corresponding to degradations of the stability by 33% and 40%, 
respectively. If the laser linewidth is larger than the natural linewidth, the probe time 
can be shortened further to adjust the Fourier transform limit to the laser linewidth. 
In the limit ri/27r Av and when r is replaced by I/Tl in the definition of S, the 
optimization yields t = 2.44/r^, p max = 0.78, and Sp — 2.71 for Rabi excitation, and 
t = 1.00/Tl, Sr = 2.59 for Ramsey excitation. The regime of a negligible natural 
linewidth is relevant if a highly forbidden reference transition like the 2 5i/2 — > 2 F 7 / 2 
octupole transition in Yb + with its nanohertz natural linewidth is excited 

The optimization problem that was discussed in this section has also been treated 
in two recent publications |12l I13| with emphasis on slightly different parameter 
regimes. Champenois et al. |12| also consider the influence of a non-zero vibrational 
quantum number of the ion and of high laser linewidth in the range Tl/2tt > Ai/n. 
For the case considered here ((n) ss and Fl/27t -c Afo), they propose excitation 
profiles (figure (5) in ^21) that differ from the ones found here (figure (2)) and yield 
instabilities that are about 50% larger. In the work of Riis and Sinclair the 
probe laser detuning that is used for the interrogation of the ion is always adjusted 
to the point where the slope of the excitation spectrum is maximal, and is not fixed 
to the half-linewidth detuning as in our calculation. For Ramsey excitation, these 
two criteria agree and identical optimum excitation parameters are found. The Allan 
deviation that is predicted in 13 on the basis of these parameters is about 35% 
lower than that calculated here using equations (7,8). This is due to the fact that 
our linear approximation Ap/ Aluq for the slope of the Ramsey excitation spectrum 
in the definition of S underestimates the maximum slope by a factor 2/tt. For Rabi 
excitation and Tl <C 1/t, the additional degree of freedom in reference p^l leads to 
an optimum with slightly shorter probe time (t = 1.63r) and higher Rabi frequency 
(ojp — 2.1PT) than found by us. The optimized Allan deviation in ^1 is about 30% 
lower than the value found here with equations (6,7). 
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Figure 3. Influence of cycle dead time on the stability for the cases of Rabi 
(solid lines) and Ramsey excitation (dot-dashed lines). With increasing dead time 
the optimized parameters shift to longer probe times (a) . The minimum stability 
parameters (b) indicate that Rabi excitation becomes more favorable than Ramsey 
excitation for > 3r. 



4. Influence of cycle dead time and the Dick effect 

A non-negligible cycle dead time reduces the rate at which the spectroscopic 
information can be gathered. Considering this parameter in the optimization, it turns 
out to be advantageous to use longer probe times in order to prevent the ratio t/t c from 
becoming too small. Figure (3) shows how the optimized probe time increases and 
how the stability degrades with increasing dead time. For the case of long dead time 
{td 3> t), the optimized probe times converge to limits of 5.07r for Rabi excitation 
and 1.99r for Ramsey excitation, while the stability parameters diverge. Since for long 
probe times Rabi excitation results in a higher excitation probability (cf. figure (la)), 
it provides a lower instability than the Ramsey method for dead times that are longer 
than 3t. From figure (3b) it can be seen that for an optimized interrogation sequence 
the dead time should be kept below w 0.2 r. In order to satisfy this requirement, the 
ion must be actively returned to the ground state after the state detection. 

A further potentially harmful influence of the dead time on the stability comes 
from the so-called Dick effect [201 12] : the long-term stability of the standard may 
be degraded by down-conversion of the frequency noise of the interrogation oscillator 
at Fourier frequencies near the harmonics of l/t c . To evaluate the influence of this 
effect the noise spectrum of the oscillator has to be known. An analytic estimate of 
the limiting instability was given for a noise spectrum dominated by flicker frequency 
noise. For the case of Ramsey excitation one finds 21 : 



where a yosc is the flicker floor instability of the oscillator. With already achieved 




(10) 
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experimental parameters like t/t c > 0.6 and a nicker floor minimum a yosc < 5 ■ 10 

it can be seen that the limitation from the Dick effect cr y lim w 2 • 10~ 16 \Jt c jT 
is well below the QPN-limited instability of a y w 10~ 15 -\/i c /T for present single-ion 
frequency standards. 



5. Servo algorithm 

A single interrogation of the ion provides only one binary unit of information that 
has to be processed by a servo system before it can be used to act back on the laser 
frequency. In this section we describe a suitable servo algorithm and investigate its 
dynamic properties. 

Suppose the laser oscillates at a frequency /, close to the center of the resonance 
line. A sequence of 2z cycles is performed in which the ion is interrogated alternately 
at the frequency / + = / + 5 m and at /_ = / — S m . The numbers of successful 
excitations n + at / + and rt_ at /_ are counted. After an averaging interval of 2z 
cycles an error signal is calculated as 

e = m , (11) 

z 

and a frequency correction g ■ e is applied to the laser frequency before the next 
averaging interval is started: 

f^f + 9-e. (12) 

The numerical factor g determines the dynamical response of the servo system and we 
will call it 'gain' for short. Since the frequency corrections calculated for subsequent 
intervals are added up, this scheme realizes an integrating servo loop. 

In order to obtain the maximum slope of the excitation spectrum, the best value 
for S m will be close to the half linewidth of the atomic resonance. The time constant 
and the stability of the servo system are determined by the choice of the parameters 
g and z. If the laser frequency / is initially one half linewidth below the atomic 
resonance and if p max = 1, the resulting value of (n+ — nJ)j z will also be close to one. 
Consequently, with g w 1, the laser frequency will be corrected in a single step. If 
g <C 1 , approximately 1 / g averaging intervals will be required to bring the frequency 
close to the atomic resonance. For g ss 2, one expects the laser frequency correction to 
jump between — S m and +S m without ever settling at a small value. If the error signal 
is averaged over a large number of cycles 2z, the response time of the servo system 
is increased and the demands on the short-term stability of the probe laser become 
more stringent. For a small value of z, however, the short-term stability of the system 
may be unneccessarily degraded by strong fluctuations in the error signal because of 
QPN. 

In order to guide the choice of suitable servo parameters, we performed numerical 
Monte Carlo simulations of the servo action under the influence of QPN. As an 
approximation of the conditions obtained with Rabi excitation, an excitation spectrum 
of the form 

, f s ( sin(7r(/- v )ha) \ 2 . . 

P{f) = Pmax 77 ri (13) 

V n/ ~ ^0)/70 / 

was used. Here, 70 determines the linewidth (70 ps 1.13Ao;o/27r) and 5 m = 7o/2 was 
chosen. The time constant of the servo was determined from the exponential response 




Figure 4. Scaled Allan deviation 0" H fo/70 of the laser frequency calculated from 
simulated realizations of the servo algorithm for different values of the servo gain g. 
The free-running laser (g = 0) is assumed to perform a random walk of frequency. 



of the error signal e(t) oc exp(— t/t servo ) to a step-like frequency perturbation and was 
found to be well described for g/z < 0.2 by 

i scrvo « 1.37 t c . (14) 

9 ' Pmax 

Since the time required to determine the error signal is 2zt c , the servo time constant 
can not be significantly shorter. With equation (14) this corresponds to a limit 
9 • Pmax < 1.46. In fact, stable operation with a time constant of about 2t c can 
be obtained for g = 1, p max = 1 and z = 1. Damped oscillatory behaviour of the 
servo is observed to set in at g ■ p max ~ 1.5. 

The Monte Carlo simulations can also be used to study the instability of the 
frequency standard that is obtained with this servo algorithm. An excitation spectrum 
like in equation (13) withp max = 1 was assumed and laser frequency noise was modeled 
as a random walk of frequency with maximal step size 7o/200 and 10 steps per cycle, 
corresponding to conditions where the short-term laser linewidth is smaller than 70, 
as realized in some experiments ^| |S] . In figure (4) the Allan deviation is plotted 
for the free-running laser (increasing as vT) and for servo parameters z — 4 and 
different values of g. For a small gain like g = 0.05, the short-term instability is 
mainly determined by the free-running laser. & y (T) then reaches a maximum at 
T w 250t c or T w 2t sorvo and decreases like 1/y/T for T > 1000i c . The value 
g = 0.05 is too small to completely eliminate the influence of laser noise at long 
averaging times, as can be seen from the comparison with the curve obtained for 
g = 0.8: Here, the short-term stability is worse than that of the free-running laser, 
but the long-term stability reaches an optimum that is only limited by QPN and 
is independent of g, z, and the laser frequency noise. This minimized instability 
follows the scaling a y — 0. 29(70/^0) ■ y/t c /T for the excitation spectrum in equation 
(13) with p max = 1. We verified numerically that for smaller values of j? max the 
numerical prefactor increases like yj (2 — p m ax) /p m ax) , in agreement with equations 
(1,7). Equation (13) can be used to approximate the shape of the spectrum for Rabi 
excitation in figure (2) with parameters 70 ~ 1/2t and p max w 0.5 for t c w 2r. 
This leads to an Allan deviation a y = 0.36/(^oV tT) that can be compared to the 
scaling obtained in section (3) for Sp — 3.22: a y — .51 / (voVtT) ■ As discussed 
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above, this indicates that the Allan deviations as calculated from Sp and equation (7) 
are somewhat pessimistic because of the approximation we used for the slope of the 
spectrum at the FWHM point. 

In a practical realization of the algorithm one has to take into account that some 
cycles may fail to provide clear information on the effect of the last probe pulse, for 
example because of momentary heating of the ion or of an extended sojourn in the 
metastable level that extends into the following cycle. We call these cycles 'invalid' 
for short. In our experimental realization, for symmetry reasons an equal number of 
interrogations at / + and /_ is processed in each averaging interval. To achieve this, 
the sign of the laser detuning is not changed after an invalid cycle, but the cycle is 
repeated at the same detuning. If at the end of an averaging interval the numbers 
of valid cycles at /+ and at /_ are unbalanced the result of the last unpaired valid 
cycle is discarded. If invalid cycles occur with a probability p- mv , the long-term Allan 
deviation is increased by a factor 1 / y/1 — p; nv because of the reduced rate of frequency 
correction information. 

In summary, the results of the Monte Carlo simulations suggest that a small value 
of cycles 2z per averaging interval should be chosen in order to optimize the temporal 
response of the servo system. Gain values in the range 0.1 < g < 1 will provide the 
best overall stability, depending on the amount and spectral characteristics of the laser 
frequency noise. 

6. Servo response to laser frequency drift 

A significant servo error may occur if the probe laser frequency is subject to drift. The 
short-term frequency stability is usually derived from a high-finesse Fabry-Perot cavity 
made from materials with small coefficients of thermal expansion like ULE or Zerodur. 
These materials show long-term dimensional changes. Due to this and as a result of 
residual thermal effects, laser frequency drift rates \df/dt\ in the range 0.05 Hz/s up to 
several Hz / s are typically observed. The actual drift rate may change significantly on 
a time scale of one hour. For a conventional first-order integrating servo not operating 
in discrete time steps, an average drift-induced error e = t seTVO df/dt is expected as the 
result of a constant linear drift [52]. Our numerical simulations yield the same result 
for the servo system considered here. It is found that the servo error is independent of 
the temporal order of the interrogations at /+ and /_ within the averaging intervals. 
Since the minimally achievable servo time constant is t seTVO « 2zt c for stable operation 
(see above), the drift-induced error can not be reduced to less than e w 2zt c df/dt. 
Thus the output of the frequency standard could differ from the resonance line center 
by several hertz under unfavorable conditions. 

An efficient reduction of this servo error is obtained with the use of a second- 
order integrating servo algorithm where a drift correction e<j r is applied to the laser 
frequency in regular time intervals <dr 



The drift correction is calculated from the integration of the error signal (equation 
(11)) over a longer time interval Td r ^> tdr 



T dr 

The long averaging time Td r is introduced to reduce the influence of QPN on the 
drift correction e^. If the two gain coefficients are related by k <g; g, the use of the 



/ 



/ + e dr . 



(15) 




(16) 
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second-order integrating servo algorithm has a negligible influence on the time constant 
^scrvo- It introduces a second, much longer time constant that is approximately 
proportional to gzt c /k. In our practical realization (see below) with t c « 0.1 s and 
\df/dt\ < 0.3 Hz/s, we use tdrift = Is, Thrift = 10 s, and k in the range 3-10 -4 -3-10 -3 . 
A numerical simulation with the values g = 0.2, k — 0.002, z — 4 and p m ax = 0.6 
showed two clearly separated servo time constants of approximately 50t c and 6500i c . 
Ideally, a second-order integrator will completely eliminate the error induced by a 
linear drift [22J - In the numerical simulations performed for the conditions given 
above, no statistically significant servo error due to a linear drift was detectable down 
to the millihcrtz level. 

7. Comparison with the experiment 

The servo algorithm described above is experimentally implemented in the 171 Yb + 
single- ion optical frequency standard at PTB 1231 1241 12*5| . The reference transition 
of this standard is the electric quadrupole transition 2 Si/ 2 (F = Q) — > 2 D 3 / 2 (F = 2) 
at 688 THz (436 nm wavelength) with a natural linewidth of 3.1 Hz. A single 
ytterbium ion is stored in a miniature Paul trap and is laser-cooled to a sub- 
millikelvin temperature by exciting the low-frequency wing of the quasi-cyclic 
hyperfine component 2 Si/ 2 (F = 1) — > 2 P\/ 2 {F = 0) of the resonance transition 
at 370 nm. The corresponding mean vibrational quantum numbers for the secular 
oscillations with frequencies of 0.7 MHz and 1.4 MHz are in the range 10-20. A weak 
sideband of the cooling radiation provides hyperfine repumping from the (F = 0) 
ground state to the 2 P 1 / 2 (F = 1) level. At the end of each cooling phase, the hyperfine 
repumping is switched off in order to prepare the ion in the (F = 0) ground state. The 
reference transition is probed by the frequency doubled radiation from a diode laser 
emitting at 871 nm. The short-term frequency stability of this laser is derived from 
a temperature-stabilized and seismically isolated high-finesse cavity made from ULE. 
In this work probe pulses with durations ranging from 30-100 ms are used, leading to 
linewidths from 10-30 Hz. From the fact that essentially Fourier-limited resonances 
with 10 Hz linewidth are obtained, we conclude that the linewidth of the probe laser is 
well below 10 Hz [S]- The excitation to the metastable 2 D 3 / 2 level is inferred from the 
resonance fluorescence level detected in the first 5 ms after reapplication of the cooling 
laser. The rate of invalid cycles (see section (5)) can be as high as 30% for a high 
excitation probability to the 2 D 3 / 2 level. The total cycle dead time needed for laser 
cooling, hyperfine pumping, and state detection after the probe pulse, is presently 
60 ms. 

To enable a quantitative study of the instability for long averaging times, we 
measure the frequency difference between two similar systems with single 171 Yb + 
ions stored in separate traps [22] ■ Both traps use the same cooling laser setup and 
synchronous timing for cooling, state preparation, and state detection within each 
cycle. The beam from the probe laser is split and two separate frequency shift and 
servo systems are employed to stabilize the probe frequencies to the reference transition 
line centers of the two ions. The frequency difference between the two systems is 
recorded once per second. For times that are larger than t ser vo the instability of this 
difference signal is independent of the temporal overlap of the probe pulses in the 
two systems. This seems to indicate that the observed long-term instability of the 
frequency difference is not affected by correlations of the servo responses to short- 
term laser frequency fluctuations. Temporal changes of the laser frequency drift rate 
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Figure 5. Allan deviation of the frequency difference between two Yb + single- 
ion optical frequency standards at 688 THz. The experimental data points are 
compared to Monte Carlo simulations (solid lines) of the QPN limited instability 
for the experimentally realized servo parameters. The reference transition 
linewidth is 30 Hz (a) and 10 Hz (b). Two solid curves are plotted in (a) because 
slightly different servo parameters were used in the two systems (see text). 



did not contribute to the observed instability, even when the time constants of the 
two second-order integrating servos were not matched. Consequently, the experiment 
allows us to measure the combined instability of two independent single-ion frequency 
standards. 

Figure (5a) shows the Allan deviation of the frequency difference for Rabi 
excitation with t = 0.03 s and t c = 0.09 s. Here the reference transition was 
resolved with a linewidth of 30 Hz in both traps, so that the influence of the probe 
laser linewidth was negligible. For averaging times T > 80 s, the experimental 
data points approximately follow the scaling o~ y (T) « 1.1 • 10~ 14 /y/T/s. Under 
the assumption of identical behaviour of the two standards, this corresponds to an 
instability cr y (T) s» 7.8 • 10~ 15 /y/T/s for the individual system, about a factor 
of three above the optimal stability that is achievable for idealized conditions (cf. 
section (3)). The solid lines show the results of Monte Carlo simulations for the 
realized servo parameters S m — 13 Hz, g — 0.2, and k = 0.0005. The shape of the 
excitation spectrum was obtained from the integration of the Bloch equations for the 
experimental pulse parameters. Two curves are presented since the two systems were 
operated under slightly different conditions: The curve with the lower Allan deviation 
is calculated with z = 5, p max — 0-65 and p lm — 0.30, while the curve with the 
higher Allan deviation represents the conditions z = 4, p max = 0.55 and p lm = 0.35. 
According to equation (14) these sets of parameters result in approximately equal servo 
time constants t S ervo ~ 7 s. Both simulated single-system instabilities were multiplied 
by y/2 for comparison with the instability of the frequency difference. 
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The Monte Carlo simulations assume a monochromatic laser excitation and 
include QPN as the only source of noise. It can be seen that the agreement between 
experiment and simulation is good and that the experimentally achieved instability 
is essentially consistent with the QPN limit, except for an excess instability of about 
20% at the longest averaging times. 

Figure (5b) shows a similar evaluation for a reference transition linewidth of 
about 10 Hz, obtained with Rabi excitation and t = 0.09 s, t c = 0.16 s. According 
to the optimization shown in figure (1), a linewidth of about 10 Hz will provide the 
lowest possible instability for the natural linewidth A^o = 3.1 Hz of the 688 THz 
reference transition of Yb + . Both systems were operated under identical conditions: 
S m = 4 Hz, g — 0.4, z = 5, p max ~ 0.35, and p lnv sw 0.3. The servo time constant 
estimated from equation (14) is 6 s. The experimental result for the instability of the 
frequency difference is (Tj,(100s) = 9 • 10~ 16 . Again, the solid line shows the result 
of the simulated single-system instability, multiplied by \[2. Here an excess noise of 
about 50% above the QPN limit is observed in the experimental data over the whole 
range of averaging times. This stability degradation might be related to probe laser 
frequency fluctuations and to uncorrelated fluctuations of the optical path lengths to 
the traps. 

8. Conclusion 

We have discussed practical aspects of the problem of the stabilization of a laser 
frequency using the signal from a single trapped ion. With an experimental 
implementation in the 171 Yb + optical frequency standard we have obtained a quantum 
projection noise limited instability below cr y (T) = 1 • 10~ 14 /-\/T/s, comparable to the 
best results obtained so far with optical frequency standards based on ensembles of 
more than 10 6 cold Ca atoms and in a comparison of standards based on a single 
Hg + ion and on a Ca atom cloud |2*7] , 

Some technical improvements of our standard appear straightforward. In 
particular, it would be beneficial for the stability to reduce the rate of invalid cycles. 
An additional hyperfine pump frequency could drive the 2 5 l 1 / 2 (-F 1 = 1) — ► Pi/ziF — 
1) component of the resonance transition during the state preparation and a quenching 
laser may be used to deplete the metastable 2 D 3 / 2 (F = 2) state after an excitation 
has been detected. These improvements may also allow us to reduce the dead time to 
about 20 ms and to approach the optimal Allan deviation of o y {T) w 2 ■ 10~ 15 / v/T/s 
that is estimated for this system. 

Quantum projection noise imposes a fundamental limitation on the stability of a 
frequency standard that uses a single ion as a reference. We have discussed how the 
interrogation parameters have to be chosen in order to obtain the lowest possible long- 
term instability, taking realistic experimental parameters into account: For a reference 
transition with natural linewidth Avq, a laser linewidth that is smaller than O.lAz^o, 
and a cycle dead time < 0.2/(27rAi/ ), an instability a y » 2{/S.v$ / vq) 1 / {2tt £\v Q T) 
will be achievable. For an optical frequency Vq rs 1 • 10 15 Hz, a very narrow linewidth 
Ai/o ~ 0.01 Hz would have to be chosen if the instability shall reach the level of 10~ 18 
within a few hours of averaging. While suitably narrow reference transitions exist in 
several ions |19l 1281 125] . substantial improvements in laser linewidth narrowing are 
required to progress into this regime of ultrahigh resolution spectroscopy. 
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